Perturbative calculation of the scaled factorial moments in second-order quark-hadron 
phase transition within the Ginzburg-Landau description 
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The scaled factorial moments Fq are studied for a second-order quark-hadron phase transition 
within the Ginzburg-Landau description. The role played by the ground state of the system under 
low temperature is emphasized. After a local shift of the order parameter the fluctuations are around 
the ground state, and a perturbative calculation for Fq can be carried out. Power scaling between 
Fq's is shown, and a universal scaling exponent v ~ 1.75 is given for the case with weak correlations 
and weak self-interactions. 
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It is well-known that ultra-relativistic heavy-ion collision is the unique way to study the vacuum properties of 
' quantum chromodynamics (QCD) in the laboratory. In the collisions the kinetic energies of the colliding particles 
CO . are converted into thermal ones, and a hot new matter state, quark-gluon plasma (QGP), might be formed. The 
system will cool with its subsequent expanding and will undergo a phase transition from the deconfined QGP to 
confined hadrons. Since only the final state particles in the collisions are observable in experiments, one may be 
asked to search for the signals about the phase transition from only those particles. Since the existence of the phase 
transition is associated with properties of the nontrivial chromodynamical vacuum, the study of quark-hadron phase 
' transition has been a hot point in both particle physics and nuclear physics for more than a decade. Besides the 
T— I ' unique features of QCD the lack of control of the temperature in the phase transition distinguishes the problem 
from the standard critical phenomena such as ferromagnetism. The nonperturbative nature of hadronization process 
in the phase transition precludes at this stage any observable hadronic predictions from first principles, and some 
. approximate models are used. One of the models is the Ginzburg-Landau model which can be used as a framework 
to calculate various moments of the multiplicity distribution and has been used in the studies of the scaled factorial 
moments in both first- and second-order phase transitions, the multiplicity difference correlators and the 

• multiplicity distributions in the phase transitions 

[ In the Ginzburg-Landau description of a second-order phase transition, the scaled factorial moments can be ex- 

• pressed as H 
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: F,^ JJfl = ^ y T^'P [J^ dz I g exp(-F[0]) (1) 

^ with Z ^ J X'0exp(-F[0]), the free energy functional F[(j>] = J dz[a \ (f) +b \ (j> \* +c \ sj(j> p], a cx (T - Tc) 
representing the distance from the critical point, b and c larger than zero. Here is associated with the hadronic 
multiplicity density of the system, and Jg dz means integration over a small bin with width 5 in the phase space. Similar 
expressions can be derived for other quantities mentioned above. In all former studies of second-order phase transition 
the gradient term in the fimctional is simply taken to be zero, i.e., the field cf) is regarded spatially uniform. 
The spatial integral of the functional over a two-dimensional bin with size 5^ is then — 5^{a | 0^ | +6 | (/) |*). 
This is of course a very crude approximation. The advantage of such an approximation is that it turns the functional 
integration into a normal one. Thus, the calculation becomes quite easy under the approximation. For a > the 
functional takes its minimum at = corresponding to the quark phase, and for a < the minimum is at > 
corresponding to the hadron phase. In all the studies the interested region is for a < 0. Numerical results do not 
show the so-called intermittency behavior, but the F-scaling, Fq oc F^' with universal scaling law Pq = [q — , is 
shown to be valid. The exponent v is called as a universal one in the sense that it is insensitive to the values of the 
parameters in the model and that it is completely determined by the structure of the functional concerned. 

The contributions from the gradient term to the moments and to the exponent v should be evaluated in some way. 
Once the gradient term is taken into the functional, one is faced with serious difficulty in the calculations, considering 
the fact that the value of parameter h for the (j>^ term can be determined in no way from first principles or from 
experimental output and may be very large. Even if the parameter h is indeed very small, negative value of a in our 
interested region also excludes the possibility of performing the usual perturbative calculations. The role played by 
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the gradient term is investigated in Q and Q . In g in each bin is still uniform, but the values of (j) in all neighboring 
bins are taken to be 4>o, field configuration corresponding to the minimum of "potential" V{(j)) = a \ cj) \^ +b \ (f) \^ . So 
the square of the gradient of is — (/)o)^- This approximation also transforms the functional integration into a 

normal one. Numerical results show that the universal scaling law f3q = {q — 1)'^ is still valid and that the exponent 
v is almost the same as without the gradient term. In the details of spatial fluctuations of in a bin is simulated 
by the Ising model for one-component spins s. Each bin is assumed large enough to contain several spin sites. This 
time, the exponent v depends on the unknown temperature, and, after averaging over the temperature, v is still in 
the range given in [|l| and |^ . 

Though the simulation in Q| is convincing, it is for lattice with one-component spins. In the Ginzburg-Landau 
model for a second-order phase transition, the field (f) is a complex number, or in other words, 4> has two components. 
At first glimps, the simulation in j^] does not correspond to the real problem discussed in the Ginzburg-Landau model, 
but as will be explained soon in this paper, it relates to the physics in an indirect way. 

In [0 , it is tempted to investigate the universality of the exponent v, with the spatial fluctuations of the phase angle 
of the complex field 4> fully taken into account. As will be shown below, the contribution from spatial fluctuations of 
the phase angle of the field (j) can be evaluated in a complete and rigorous way, and the integration over the spatial 
fluctuations of the phase angle of the field cf) will reduce the problem to one with a one-component field. 

The first observation is that all terms except the gradient one in the functional integral of Eq. (|l|) depend only on 
I (/) p. Then it is convenient to write the two-component field as a complex number in the form (j> = 0R exp(«0i). 
The spatial fluctuations of the field can be those of the magnitude (jjR and/or of the phase angle 0i (or orientation in 
an abstract space). The gradient term turns out to be 

I V0 l'= (V0r)' + -^rXv^i)' ■ (2) 

Generally, the phase angle 0i can be in any form, and the full contribution from its fiuctuations must be evaluated. 
Fortunately, the integral over can be carried out easily since it is of Gaussian form. Then one transforms the 
two-fold functional integral into a one-fold one and gets 

with functional exactly the same form as the original F[(j>]. The important difference between this expression 

from Eq. (|^) is that the functional integral variable in this new expression is a real function instead of a complex 
function as in Eq. (|l|). Then fq and Fq can be simulated by a one-component field as in Ref. 

Now we take the field (/jr (magnitude of 0) uniform, or in other words, the gradient term of 0r is omitted. 
(Galculations based on this approximation will be referred to mode 2 in this paper.) Based on the work Ref. Q one 
can drop off the V0R term, because the problem now is exactly within the scope of Ref. ||^, and the conclusions in 
Ref. 1^ encourage us to neglect the spatial fluctuations of (/)r as long as the universal scaling exponent is concerned. 
Then one gets the factorial moments as functions of variable x 



' Jo" dy expixy^ ~ y'^) 



with X = aS^/^ /b^^'^. From this expression the scaled factorial moments InF^ can be calculated, and the results are 
shown as functions of — In a; in Fig. ^ for q from 2 to 8 within the range x G (0.5, 4.0). One can see clearly that no 
strict intermittency can be claimed since all Fq approach finite values in the small x limit. So, no intermittency is 
shown in the phase transition, as shown in former studies. More importantly, the power law can be found between 
Fq and F2, as shown in Fig. ^ with the same data as in Fig. |l|. 

For the convenience of comparison with former case, we write down the expressions of the scaled factorial moments 
without spatial fiuctuations (mode 1 in this paper), which can be read 

J, _ fi^ dyy"} expjxy - y^) 
dy exp{xy - y'^) 

with X = aS/Vb. Numerical results for InFq in this mode are shown in Fig. |^. In the upper part of the figure \nFq 
are shown as functions of — In a; for q from 2 to 8 with x in the same interval x € (0.5,4.0), and in the lower part 
InFq are shown as functions of I11F2 with the same data as in upper part. One can see from upper part of the figure 
that the general behaviors of In Fq as functions of — In a; is similar to those in Fig. |^, though the definition of x in 
this case is different from that for Fig. |l|. The values of Ini^^ in the two cases are also different. For same value of x, 
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InFg in the former case have larger values. This difference is reasonable if one notices the difference in the definition 
of variable x. What interests us is the scaling law between Fq and F2. The power law scaling between Fg and F2 can 
be seen obviously in the lower part of Fig. pi the same as shown in other studies cited in the references. 

From Fig. |2| and the lower part of Fig. one can get the scaling exponents l3q for the two different modes by 
fitting the curves. Pq can also be given analytically. One can expand the expressions for In Fq in the two modes as 
power series of x in small x limit, and then one gets the slopes Kq for In Fq and /3q = Kq/K2- The expressions for Kq 
for the two modes in this paper are 



r 0/2 + 1) r(3/2) , , r i 

Kq= ' -q ' +{q-l) , ' for mode 1 

' r(g/2 + i/2) ^ r(i) r(i/2) 

r(g/2 + 3/4) £(5/4) . r(3/4) 

r(g/2 + i/4) '^r(3/4)'^^^ ^r(i/4) 



- T - + (g - l)^, formode2 



One can find only a small difference between the exponents v from these two expressions. The results are shown in 
Fig. ^. In mode 1 (without spatial fluctuations) 1.3335, and in mode 2 (with spatial fluctuations of the phase angle 
of the field (j)) v=1.2n2. The exponents obtained from these analytical expressions are very close to the ones from the 
fitting. The universal exponents v are also very close to one another and can be regarded as the same within accuracy 
4%. Physically, these two modes correspond to different situations. In mode 1 no spatial fluctuation of (p is in the 
problem, but in mode 2 the spatial fluctuations of the phase angle of the complex field (p are fully evaluated. Since 
these two different considerations give very close exponents one can say that the exponent v is indeed insensitive 
to the spatial fluctuations of the phase angle. 

It is, of course, very interesting to investigate directly the effect of the term (v0r)^ on. the moments, which is the 
main topic in this paper. 

Our second observation is that the final state particles are in a finite phase space at any high but finite colliding 
energy. This means that the fluctuations of the field (p should not be uniform since the field must be zero in the 
region excluded by the conservation laws. Thus there exists a boundary condition for 0R. For the convenience, we 
use (j) instead of (/ir in the following of the paper if no confusion will be aroused. The boundary condition of cj) is of 
Dirichlet type in our problem because of the fact that the particle density out of a finite region should be zero. In 
the following, we only discuss a one-dimensional phase space such as the rapidity, and the boundary condition can, 
not losing any generality, be written as 0(0) — 4>{L) = 0, with L the length of the finite phase space interval. With 
the gradient term in the functional, the functional integral can only be calculated perturbatively. But there are two 
important differences from the usual perturbations. The first difference is the finite size of the phase space. The 
second is the non-positivity of the coefficient of the Gaussian term in the functional F[(j>]. So, some new techniques 
are needed which will be discussed in this paper. 

The organization of the paper is as follows. In Sec. II we discuss the ground state of a finite-size system under various 
boundary conditions. In Sec. Ill a new perturbative calculation scheme is proposed with the effect of local spontaneous 
symmetry breaking taken into account. In Sec. IV we calculate the scaled factorial moments perturbatively. Sec. V 
is for our main results and conclusions. 



II. LOCAL SPONTANEOUS SYMMETRY BREAKING FOR FINITE-SIZE SYSTEM 



Finite-size effects near critical points have been remained over the past two decades to be an important topic of the 
active research both theoretically and experimentally ||^ in condensed matter physics. Nowadays, the experimental 
sample are usually so pure and so well shielded from perturbing fields that the correlation length can grow up to 
several thousand angstroms as the critical point is approached. When one or more dimensions of a bulk system is 
reduced to near or below a certain characteristic length scale, the associated properties are modified reflecting the 
lower dimensionality. It is believed that flnite-size effects are precursors of the critical behavior of the inflnite system 
and can be exploited to extract the limiting behavior. A central role plays the finite-size scaling behavior predicted by 
both the phenomenological § and rcnormalization group ||l^ theories. Those theories allowed a systematic discussion 
of the finite-size effects and, consequently, form the cornerstone of our current understanding of the way in which the 
singularities of an infinite system are modified by the finiteness of the system in some or all of the dimensions. Of 
course, the exact form of scaling functions can't be given in those scaling theories. 

In 1985, Brezin and Zinn-Justin (BZ) |0 and Rudnick, Guo and Jasnow (RGJ) ||l2| developed two field-theoretical 
perturbation theories for the calculation of the finite-size scaling functions within the 0^ model which corresponds to 
the Ising model. Most applications of these theories to three-dimensional systems have been restricted to T higher 
than the bulk critical temperature Tc [ p^ with a few calculations in region below Tq p^ ]. However, some limitations 
exist in the theories of and [Q. As pointed out in the first paper in [|l^, the theory of BZ is not applicable for 
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T < Tc and the results from RGJ theory are not quantitatively reliable in the same ternperature region since the 
coefScients of the Gaussian terms in the integrals are negative for those temperatures. In [[l6| the order parameter is 
expanded into a sum of eigenfunctions of for various boundary conditions. Again, the functional integral is turned 
out into a product of normal integrals. But the fluctuations can be evaluated only for temperature not too far below 
the critical point. Authors of tried to avoid the difficulty mathematically, but they failed to account for the origin 
of the difficulty physically. Although the modified perturbation method in can be used for both T > Tc and 
T < Tc, the calculation is lengthy and can be done only to first order in practice. Since one does not know the exact 
order of values of higher order terms, theoretical results may have large uncertainty. 

It has not been answered that which physical effect causes the failure of direct perturbative calculations of fluc- 
tuations for finite-size systems with temperature below Tc- In our opinion, the real origin of the difficulty lies in 
the lack of knowledge about the spontaneously symmetry breaking for finite-size systems. It is well-known that an 
infinite system will have non-zero mean order parameter ipo, which is called ground state of the system in this paper 
since it corresponds to minimum of the Hamiltonian H, if the temperature is below the critical one, and everyone 
knows that the difficulty of negative coefficient of the Gaussian term can be overcome by shifting the order parameter, 
— » (j) + (j)o- This phenomenon is known as the spontaneous symmetry breaking because of the fact that 0o does not 
have the same symmetry as H does. This kind of spontaneous symmetry breaking for infinite system can be called 
global since the shift 4>q is the same constant for every point in the space. For a finite-size system, such a simple 
shift of the order parameter does not work because of the existence of specific boundary conditions for the system. 
Anyway, fluctuations of the system, in their own sense, should be around certain ground state which corresponds to 
the minimum of the Hamiltonian H, and they can be approximated by Gaussian terms in most cases if they are not 
very large. Thus one sees that the ground state plays an determinative role in the study of fluctuations in the phase 
transitions at low temperature. For inflnite system, the ground state (jjg is constant and can easily be calculated. But 
for a finite-size system, the ground state is usually not a constant but depends on the boundary conditions imposed 
on it. This is understandable. For infinite system the ground state is determined completely by the self-interactions 
of the field. In other words, the ground state is dictated only by the "potential", and there is no boundary effect. 
For a finite-size system, however, the effect of the boundary must be taken into account. For the case with local 
interactions, the effect is realized through the gradient term. Thus the ground state for system with finite size is 
determined by the gradient term and the "potential" . Then the shift of the field at a point depends on the position 
in the space. So, the spontaneous symmetry breaking for finite-size system can be called a local one. Therefore, the 
solution for the ground state is non-trivial but necessary, and one has reason to hope that the difficulty mentioned 
above for finite-size systems can be overcome once the ground state is known. 

It should be pointed out that all perturbation theories mentioned above are based on Fourier decomposition of the 
order parameter. This method is natural because the decomposition enables one to transform the functional integral 
into an infinite product of tractable normal integrals. Although such a decomposition has simple physical explanation 
which is very fruitful for the understanding of properties of infinite systems and can deduce reliable physical results, as 
in the case of usual field theories in particle physics, it brings about a great deal of calculations for finite-size systems. 
This is not surprising. As is well-known, quantities complicated in coordinate space may have simple momentum 
spectra thus look simple in momentum space, but those obviously nonzero only in a finite range must have puzzling 
momentum spectra. Therefore, for the study of properties of finite-size systems, calculations in coordinate space 
might be simpler and more effective. The point here is that one must calculate the complicated functional integral 
which is very difficult to be evaluated directly. 

In this section, we first calculate the ground states for a (j)* model of a second-order phase transition with one- 
component order parameter under various boundary conditions. All the boundary conditions are useful in the study 
of condensed matter physics. Then, with the ground states, the Hamiltonian of the system is reexpressed as Gaussian 
terms and higher order perturbations of a locally shifted order parameter. And it is shown that the perturbative 
calculation can be done with the new Hamiltonian for temperatures far below the bulk critical point. 

In the (j>* model for a second-order phase transition in condensed matter physics with a one-component order 
parameter, the partition function can be expressed as a functional integral of exponential of the Hamiltonian H of 
the system 



2^ 



u 

4!' 



(6) 



in which 7 ~ a'{T — Tc), a' and u are temperature dependent positive constants, (f> is the order-parameter of the 
system. In the following, we are limited only to a film system with thickness L. Since we are interested only in 
the temperature region T < Tc or 7 < 0, the Hamiltonian H can be standardized by introducing correlation length 
C = \/— 1/7, new order-parameter ^' = (j)/4'o, with 0o = ^Z— 67/u the vacuum expectation of the order parameter for 
bulk system, scaled coordinates r' = r/L, and reduced thickness I — L/£^, into 
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H = / d^r' 



From this expression one can get the equation for the ground state by ^ = 0. The ground state ^'o(^) satisfies 



(7) 



f d'^o 



= -*o + *o ■ (8) 



fn the equation we have used z instead of z' in the range (0, 1) to denote the coordinate along the thickness direction. 
Derivatives in other directions do not appear in the equation since any state with non-zero derivatives in other 
directions does not correspond to the minimum of H . But if the system is fully limited in all directions, last equation 
should have in place of d^/dz'^. In last equation is solved analytically for Dirichlet boundary conditions 
^-(0) = *(1) = 0. The exact solution is 

*o(z) = -^^sn(2zF(fc),A:) , (9) 
V 1 + k"' 



in which k is determined by / through I ~ 2vT+fc^F(fc). Here, F{k) is the first kind complete elliptic integral, sn(2;, k) 
is elliptic sine function. Unfortunately, no simple compact solution is found yet for other boundary conditions. One 
can easily see that the main obstacle comes from the nonlinear term in the second-order differential equation of 
^'o in Eq. (H). To find approximate solutions of ^'o for other boundary conditions, the following method can be used. 
First of all, we replace in Eq. ^ by A^'o and get a solution satisfying the same boundary condition. For Dirichlet 
boundary conditions, the solution is 

'J'o = ^siuTTX , with A = 1.0-7rV^^- (10) 

The constant A can be determined by requiring the mean square of the deviation caused by the replacement, i.e., the 
integral dz{^Q — A^'o)^, to be minimum. Thus one gets 



siuTTZ . (11) 



Now one can see that the requirement of a minimum deviation caused by the replacement is equivalent to retaining 
siuTTz term but neglecting terms with higher frequency in VPg. Thus, this approximation is equivalent to the standard 
functional variation method. The virtue of this method is that it can be used simplier and in a step-by-step way. As 
discussed in ||l^ the ground state is ^'o = if the reduced thickness I of the film is less than n. The existence of 
minimum reduced thickness of the film implies a shift of the critical temperature for the finite system from the bulk 
one. The exact solutions and the approximate ones are compared in Fig. ^ for I/tt—I.OB, 1.10, 1.15, and 1.20. A 
very good approximation can be seen. For larger /, the same approximative method can be used further after shift 
= + i/4(l - TT^/P)/ismTTz in Eq. (||). For Neumann boundary conditions, *o(0) = ^'o(l) = 0, the ground 
state can also be obtained in a similar way. The result is 

^-0(2) = 1.0 for T<Tc (12) 

Then one can consider mixed boundary conditions ^'o(O) = 0,^*9(1) = 0. The first order approximation of the 
solution for ground state is 




*o(^) = W3 (^1-^jsin— for l>n/2. (13) 
As a final example, we give the ground state for periodic boundary condition '^oiz) = '^o{l + z). The ground state 



«'o(z) = 1.0 for T<Tc (14) 

Though the ground state for periodic and Neumann boundary conditions are the same the fluctuations of the fields 
in the two case are different. It should be pointed out that — V&o is also a ground state of the system. Then the 
fluctuations of the system can be around either '$0 or — ^o- This is the copy for finite-size systems of spontaneous 
symmetry breaking in cj)'^ model. The significant difference from the usual spontaneous symmetry breaking is that the 
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ground state is usually not a constant and depends on the boundary conditions. So that we have a local spontaneous 
symmetry breaking in this paper. With the ground state ^'o, one can locally shift the order parameter 5* = 5*' + 4*0, 
then the Hamiltonian H turns out to be 



e 



(15) 



In this expression, i/['I'o] has the same form as H[$] in Eq. with ^'o in place of 'J. Now the quadratic part of 
fluctuation ^f' is positive definite for I larger than a characteristic length, or for temperature enough below the critical 
point. Then one sees that the new Hamiltonian can be safely used to calculate perturbatively fluctuations at low 
temperature region for finite-size systems. Then the difficulty of the negative coefficients of the Gaussian terms is 
avoided after the effects of local spontaneous symmetry breaking are taken into consideration. 



III. PERTURBATIVE THEORY FOR FINITE-SIZE SYSTEM UNDER T < Tc 

FromEq. (|l|), a new perturbative theory can be developed for finite-size system with local spontaneous symmetry 
breaking. First of all, one can introduce for an one dimensional system a generating functional Z[J] 



Z[J] = j X>0exp(-i/ + j dzJ(f>) 



(16) 



The generalization to more general cases is obvious. Up to an unimportant constant factor, the generating functional 
for a one-dimensional system can, in a standard way, be written as 



Z[J] = exp(Ai J dzJ"^Q)eyij>i^~X J dz 



XiSJ 



XidJ 



exp 



IM J dzdyJ{z)G{z,v)J{y) 



, (17) 



with Ai = LyJ&yf\/u ~ L(f>Q, X ~ %L^I^ ju. In last equation, the Green's function G'(z, y) satisfies 



1^ 

^2 



l + 3*^(^) 



G{z,y)^8{z-y) 



(18) 



The first factor in the generating functional shows a great difference between present theory and the usual ones that 
there exists a nontrivial solution for the classical equation 5H / 5(j) = J for J = 0. For systems with higher dimension 
d> 1 the only changes are with L'^ in place of L in the expressions for parameters A and Ai and with place of 

(P jdz^ in last equation. The Green's function G{z,y) describes fluctuations in the full space and determines how the 
fluctuations at different points are correlated. If one can get the solution for G{z,y) for higher dimensional system, 
the fluctuations can be evaluated in the same way as for one dimensional system. Thus in the following we do not 
distinguish one and higher dimensional systems, and dz is used to represent the integral element over a volume in 
certain space. From Eq. (p^, it can be seen that each Green's function G is associated with a factor 1/A. Ai can be 
regarded as a factor associated with the external source field J. Since the derivative terms in the second factor in Eq. 
( p7| ) with respect to the external source field J will generate terms with more factors of G in the generating functional, 
the contribution of them is small if the parameter A is big enough. Then those terms in the generating functional can 
be regarded as perturbations. From the expression of A it is clear that a large A is equivalent to a small u for fixed 
L and 7. Thus the condition of a large A is consistent with that in usual perturbation theory. Then one has all the 
four ingredients diagrammatically represented in Fig. || for the perturbative calculations with the Feynman rules: 

(a) the ground state: Xi'^o(z), 

(b) the Green's function (propagator): -YG{z,y), 

(c) three-line vertex: — / dz'^ oi^), 



(d) four-line vertex: 



^^Jdz. 
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Using these ingredients all physical quantities can be calculated. For example, to the first order of the perturbations, 
one has 



= *o(^) - I y" du'ifo{u)G{u, u)G{u, z) , 

{{^{z) - Mz)) {^{y) - My))) = \G{z,y) j duG{z,u)G{u,u)G{u,y) . 

Here, the symbol (• • •) represents the average over the fluctuations, the range of the integral over u is (0, 1). 

A most important feature of the perturbation theory is that all the calculations can be done in coordinate space. 
Once the non-trivial ground state is known, one can get the Green's function (propagator) G{x,y) from Eq. jl^), 
and other quantities can be obtained from Eq. ( p^ by directly taking derivatives with respect to the external source 
field J. This scheme can be used in calculating properties of finite-size systems in condensed matter physics for 
temperatures T <C Tc. 

Next section will calculate the scaled factorial moments in a second-order quark-hadron phase transition as an 
example of the applications of the perturbation theory. 

IV. THE SCALED FACTORIAL MOMENTS IN THE GINZBURG-LANDAU MODEL 

Now we turn to the calculation of the scaled factorial moments in Eq. (|l|) in a second-order quark-hadron phase 
transition within the Ginzburg-Landau description. In this description, the free energy functional F[(j)\ is in place of 
the Hamiltonian H in last two sections. After integrating over the phase angle of the field the functional remains the 
same form with a real 0r in place of the complex as discussed before. Although there are very important differences 
between normal phase transitions in condensed matter physics and a quark-hadron one, the mathematical form in the 
Ginzburg-Landau description for them is the same. In the Ginzburg-Landau description for a quark-hadron phase 
transition, the integral variable z is not in coordinate space but represents a collection of measurable quantities such 
as rapidity and azimuthal angle etc. In the following, z is identified to the rapidity. For such an one dimensional 
system, the local spontaneous symmetry breaking is also given as in Sec. II. A generating functional can also be 
introduced in the same way as in last section. The only changes are the expressions for the parameter A, Ai, and 
Here we only mention the expression for I. In present case, the correlation length is ^ = y^c/H, so I = L^\a\lc. The 
parameter c has a simple physical meaning. From the free energy functional one sees that the correlation between 
fields at different points is realized by means of the gradient term. If c is small there is weak correlation between the 
fields at different points. Thus the effective length I can be used to measure the strength of the correlations for fixed 
L and \a\. For a system at fixed temperature c is small if there is weak correlation, and vice versa. When Z — > oo, 
one may expect that the influence of correlation can be neglected and that the effect of boundary condition can be 
neglected. In the calculation of the scaled factorial moments, the factor Ai will be cancelled. So Ai can be taken to 
be 1.0 in present calculations. For any parameter I the scaled factorial moments can be rewritten from Eq. (1) as 



In this expression, J^- dz represents an integral over a range of length 8. In our calculation, the integral range is chosen 
around the center of the interval (0, 1), or in other words, in the range (l/2-i5/2,l/2-|-(5/2). As discussed in the second 
last paragraph in Sec. I, the boundary condition for our case is of Dirichlet type. So the ground state ^'o is given by 
Eq. and G{z,y) is calculated from Eq. ([T8|). 

A. Zero order approximation for fq 

We first calculate the scaled factorial moments Fq in a second-order quark-hadron phase transition at the zero 
order (or tree-level) approximation to the functional Eq. (|l^). At this level the second factor in Eq. (|l^) gives a 
factor 1.0. In the expressions of/g there are contributions from g-particle correlations represented diagrammatically 
by connected diagrams in Fig. |7| and the contributions from fewer particle correlations which can be represented by 
products of disconnected diagrams. We denote the contributions to fq from connected diagrams which give the 
contribution from the pure g-particle correlations to fq. Then the factorial moments fq at tree- level can be written as 
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f tree 
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ftrcc 
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(/f + 10(/f)3/| + 10(/f)2/| + 5/f/l + lO/l/l + /, 



(20) 



(/f + 15(/f)VI + 20(/f + 15(/f)2/4- + 6/f/5= + 10(/|)2 + 15(/|)3 + 15/1/1 + 60/f/l/l + /, 
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For the connected contributions to /g'°°, there are only two topologicaUy different diagrams, as shown in Fig. |^. For 
the first type diagram with two crosses representing the ground state, the number of identical terms is A^^ = 2''~^q\. 
The factor q\ comes from the exchange symmetry of the q particles, 2'^ from the two-lines from each point representing 
a particle, and a factor 1/2 from the identities of terms with reversal order of the g-points. For the second type of 
diagrams with no cross, the number is = N^/q = 2^'^~^\q — 1)!. To calculate the diagrams, it would be useful to 
define 



dxi dx2 



dx, G{z,xi)G{xi,X2) ■ ■ ■ G{x,,y) , 



which satisfies a recursive relation 

9iiz,y) 

Is JS 
Then the contribution from each connected diagram for fq can be written as 



dugi^i{z,u)G{u,y) ^ / dzG{z,u)g^^i{u,y) 



(21) 



(22) 



first disgram: 
second disgram: 
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dz dy^a{z)gq-2{z,y)'^o{y) , 



dzgq^i{z,z) . 



So that 



^ / dzdy ^'o(2;)gg-2(^;,2;)*o(y) + ^ / dzgq^i{z,z) 



(23) 



B. First order approximation for /, 



Now we discuss fq at the first order (1-loop level) approximation of the second- factor in the functional of Eq. (|1 
At this approximation, the factor from the second term of the equation is 



Zi[J] = 
in which {GJ)^ 



1 



dz 



J duG{z 



*o(^) 
u)J{u) 



3G{z,z){GJ), 



{GJ)l 



■iG'^{z,z) + ^G{z,z){GJ)l + 
A 



{GJ)\ 



A2 



Z[J] = Zi[J] exp( / dzJ-^o) exp 



dzdyJ{z)G{z, y)J{y) 



From the functional at this approximation 
_2A , 

the factorial moments fq can be directly calculated by using Eq. ([l9|). There are many terms contributing to fq, 
among which the most interesting terms are those represented by connected diagrams in Fig. ^ with one bulb which 
is the vertex for the perturbative interactions. The sum of the contributions from the diagrams to fq will be denoted 
by /^°°P in this paper. Then up to the first order approximation of the generating functional, the factorial moments 
fq are 
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loop 
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j-tioc floop . rloop 
Jl Is + Ji 1 

f io/*'-^v3°°p + 5/r<="/i°°p + /5°°p , 

f 20/*'-^V3°°P + IS/s*'""/!™" + 6/f"V5°°P + /6°°P 



(24) 
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For the perturbative calculation to have high accuracy, we choose the parameter A large enough to guarantee the 
following two conditions: (1) / dz\{'i/{z)) — 'iio{z)\, the integral of absolute deviation of the mean value of the order 
parameter 5* from is not larger than 0.05; (2) \Z[0] — 1| is no more than 0.05. These two conditions ensure the 
contributions from higher order terms from the second-factor in Eq. can be safely neglected. So our calculations 
arc limited to only the first order approximation. Of course, higher order approximation can be made without difhculty 
in principle, only with more diagrams drawn and evaluated. In numerical calculation, S is chosen for — In 5 in the 
range (1,4). 



V. MAIN RESULTS AND DISCUSSIONS 



As discussed in last section, we choose the parameter A to be a large number to ensure the small influence of the 
perturbations. A large A corresponds to a small correlation function G{x,y) / X. So in the Ginzburg-Landau model 
for a second-order phase transition under some choice of the parameters, there is weak correlation between the fields 
at different points together with weak self- interactions. Due to the choice of a large A (Ai = 1.0) the ground state 
^'o (whose square is the hadronic density at the state) will play a dominant role in Eq. ( |l7| ) for I large enough, and 
the Gaussian and higher order fluctuations can only bring about some small corrections to the generating functional. 
Then with the choice of A we are dealing with a case with small fluctuations. Because of the large value of A the 
first term in the brackets of Eq. (|2^) plays an important role if the ground state is obviously nonzero for larger 
parameter /. Since the powers before the brackets of Eq. ( p3| ) will be cancelled, the order of the ratios fq/{fiY is 
A-f'-i), thus very small. fl°°^/{fi)'^ have the order A even smaller. Then the scaled factorial moments Fq are 
very close to 1.0. This expectation is confirmed in numerical calculations. Numerical results show that hiFq, though 
very small, have quite complicated behaviors. They increase for — In 5 within (1.5, 2.5) and then decrease with the 
increase of — In 5, as shown in Fig. ||for parameter / = 2.637r. Thus there is no intermittency in the phase transition. 
For other choices of Z/tt 3> 1 similar results can be obtained. A more important and more interesting phenomenon is 

the power scaling between Fq and F2, Fq oc F^'' , which can be expected from the similar behaviors of In Fq in Fig. 
^ and are shown in Fig. |l^ with the same data as in Fig. |. (3q can be obtained easily from a linear fitting to the 
curves in Fig. |l^. As in former studies of Fq in Refs. [|l],|j in the phase transitions, f3q satisfies a universal scaling law 

Pq = iq- ly , with v = 1.7539 for l/n = 2.63 , (25) 

as shown in Fig. |l^. In this case the universal exponent v depends only on the value of parameter I which is a 
function of parameters \a\ for the temperature and c for the correlation strength. The dependence of the exponent v 
on temperature is consistent with Ref . |^ . But the exponent v is very different from those exponents given in former 
studies. The discrepancy is caused from the different assumptions made in former studies and in present one. In 
former studies, the effect of the gradient term is neglected, but the term (which describes the self-interaction) is 
emphasized. In those studies, the factorial moments fq can be written as 

/•OO 

dyy"^ exp(xy - y^)/ / dy exp{xy - y^) , (26) 



Jo 

in which a; is a parameter representing the bin width. From this expression one can discover that the 0^ term, 
corresponding to the —y^ term in the exponentials, is very important and cannot be treated as perturbation for any 
parameter x. It is the term that makes the integrals finite. In present calculations, the role played by the 0"* term 
is much less important. Its function is to provide a nontrivial ground state around which the fluctuations are. 
Then that term is treated as a small perturbation and is very weak indeed with our choice of parameter A. In former 
studies the fluctuations are arouns 0o — 0. Then the discrepancy between present study and former ones can be 
understood because they belong to different physical regimes. Former studies are in the nonperturbative regime with 
trivial ground state, but present study in the perturbative one with a non-trivial ground state. 

The dependence of the universal exponent v on the parameter / is also studied for l/n > 1 in which there exists 
a non-trivial ground state. The result is shown in Fig. |2[ For long correlation length (l/n a little larger than 1.0) 
the fluctuations in neighboring bins are correlated. For these / the values of are also small, so the two terms in 
the bracket in Eq. (n9^ may have comparable contributions to fq. In this region v is quite large (about 2). With 
the increase of / the correlation between the fluctuations in neighboring bins becomes weaker and weaker, and the 
exponent ly decreases first rapidly and then slowly. When I/tt > 2.5 v approaches a constant, about 1.75. The constant 
can be anticipated by considering a case with the very weak correlations among particles more than 2 (considering the 
factor 1/A accompanied with the Green's function G{x, y)). Then if only the effects of a weak two-particle correlation 



9 



is considered, one has 
and then 

F, = i.o+cifi/{ri)' . 

Here are the binormial coefficients. Since the ratio fi/ifi)^ is assumed very small one gets InFg ~ Cgfi/ifi)"^' 
so the linear relation between Ini^^'s can be verified, and one can get 

Pg = C'g=qiq-l)/2. 

From this expression, one gets the exponent f = 1.7550, very close to the one obtained in this paper. 

As a summary, the spatial correlation of the fluctuations in a second-order quark-hadron phase transition is con- 
sidered in this paper within the Ginzburg-Landau description. We deal with a case with finite phase space and with 
negative coefficient of the Gaussian term in the functional. Because of the finite size of the space, calculations in usual 
space are simpler and more effective. Due to the negative coefficient of the Gaussian term in the functional a local 
spontaneous symmetry breaking (or non-trivial ground state) exists for finite size system. Wc emphasize on the impor- 
tance of the ground state of the system, which is a version of spontaneous symmetry breaking for finite-size systems. 
Then a new perturbation scheme is developed which is expected to be applicable in the low temperature region in the 
(/)^ model for second-order phase transitions in condensed matter physics. Then as an example of the applications, 
the scaled factorial moments Fg in a second-order quark-hadron phase transition are calculated perturbatively. Power 
scaling laws between Fg^s are shown and a universal exponent v is given. 
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Figure Captions 



FIG. 1. Dependences of InFq on the bin width — In a; after the contribution from spatial fluctuations of the phase angle of 
the f ield fully taken into account (mode 2). Curves from lower to upper are for q from 2 to 8, respectively. 



FIG. 2. Scaling behaviors of Ini^, vs lni<2 for the same data as in Fig.l. 

FIG. 3. Upper part: InFq as functions of — Inx without spatial fluctuations (mode 1); Lower part: Scaling behaviors 
between InFq and ln_F2, with the same data as in the upper part. 



FIG. 4. Scaling behaviors of In f3q as function of ln(g — 1) for the two modes. 

FIG. 5. Comparison between exact solutions and approximate ones for Eq. ^ under Dirichlet boundary conditions for 
//7r=1.05, 1.10, 1.15, and 1.20. The solid curves correspond to exact solutions, dotted curves are drawn according to Eq. (^l|). 

FIG. 6. Feynman diagram representations for (a) the ground state, (b) the propagator, (c) three-line vertex, and (d) four-line 
vertex. 

FIG. 7. Connected zero-order diagrams for the contributionsr to fq. In the diagrams the number of dots is equal to q and 
an integral over the coordinate in a range with length S is implied. So (a) and (b) are for /i, and (c) and (d) are for fq with q 
dots in the diagrams. 



FIG. 8. Connected first-order diagrams for the contributions to fq for q=l, 2, 3, 4, 5, 6, respectively. 

FIG. 9. Dependences of the scaled factorial moments InF^ on the bin width — In 5 from 1.0 to 4.0 for parameter I = 2.637r 
for q—2, 3, 4, 5, and 6. 



FIG. 10. Power scaling between Fq's with the same data as in Fig. 



FIG. 11. Universcal scaling between and {q — 1) for parameter I — 2.637r. 



FIG. 12. Dependence of the universal exponent v on parameter Z/vr. 
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